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I.  INTRODUCTION 


The  feasibility  of  modeling  the  response  of  a generic  structure 
excited  by  electron  beam  energy  deposition  and  of  correlating  the  results 
with  available  test  data  was  undertaken  by  the  Ballistic  Research 
Laboratory  at  the  request  of  the  Harry  Diamond  Laboratories.  Success 
in  this  endeavor  would  enhance  the  basis  for  utilization  of  a large- 
scale  electron  beam  facility  for  simulation  of  the  effects  of  intense 
x-rays  on  missile  and  re-entry  vehicle  structures,  thus  reducing  the 
requirement  for  costly  and  time-consuming  underground  nuclear  tests. 

l-xperimental  data  were  obtained  for  aluminum  alloy  cylinders  of 
various  wall  thicknesses  by  HDL  personnel  using  the  AURORA  facility  in 
the  electron  beam  mode.  These  cylinders  had  thin  lead  overlays  selected 
so  as  to  absorb  sufficient  energy  from  the  electron  beam  to  be  partially 
melted  or  vaporized  and  to  impart  an  impulse  to  the  aluminum  cylinders 
as  they  were  blown  away.  Films  made  with  a high  speed  framing  camera 
showed  that  this  technique  was  successful.  Large  deformations  were 
induced  in  some  of  the  aluminum  cylinders.  Many  of  these  structural 
cylinders  were  instrumented  on  the  inside  with  strain  gages  to  record 
transient  strains  at  various  positions  around  the  circumference. 

The  modeling  aspect  of  this  investigation  has  similarities  to  a 

i o 

previously  reported  BRL  study  1 , the  principal  differences  being  that 
the  earlier  work  was  concerned  with  x-ray,  rather  than  electron  beam, 
deposition  and  the  previous  structural  configuration  was  an  initially 
flat  plate.  In  the  following  sections  of  this  report  necessary  details 
of  the  AURORA  tests  will  be  provided,  the  procedures  followed  in 
modeling  this  experiment  will  be  discussed,  and  the  correlation  between 
predictions  and  test  data  will  be  evaluated. 


II.  THE  EXPERIMENTAL  CONFIGURATION 

The  structural  cylinders  each  had  an  outside  diameter  of  152.4  mm 
(6.0  inches)  and  were  fabricated  from  6061-T6  aluminum  alloy  by 
rolling  and  butt  welding.  The  stress-strain  diagram  obtained  by 
pulling  a tensile  specimen  of  this  material  is  shown  in  Figure  1.  These 
cylinders  were  190.5  mm  (7.5  inches)  long  and  were  supported  at  the  ends 


1N.  J.  Huffing  ton,  Jr.  and  H.  L.  Wisniewski,  "Structural  Response 
Induced  by  X-Rays,"  Proceedings  of  the  Fifth  Symposium  on  Nuclear 
Survivability  of  Propulsion  and  Ordnance  Systems,  October  1975. 

2 H . L.  Wisniewski,  "Coupling  of  X-Ray  Deposition  to  Structural  Response," 
U.S.  Army  Ballistic  Research  Laboratories  Memorandum  Report  No.  2761, 
June  1977.  (AD  #A042667) 
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used  for  Cylinder  Specimen 


by  stepped  discs  as  shown  in  Figure  2.  In  turn,  these  discs  are 
presumed  to  be  rigidity  secured  to  an  immovable  foundation  during  the 
test.  The  free  length  of  the  cylinders  between  the  end  discs  was 
152.4  mm  (6.0  inches);  however,  since  the  same  end  fittings  were  used 
regardless  of  the  gauge  of  the  cylinders,  it  is  presumed  that  there  was 
some  "slop"  between  cylinder  and  end  discs  for  the  thinner  gauge 
cylinders.  This  poses  a problem  &s  to  the  proper  modeling  of  boundary 
conditions  for  the  cylinders  in  the  response  analysis. 

Rather  than  consider  a spectrum  of  test  data,  the  succeeding 
discussion  will  consider  only  a single  representative  test  (Pulse  1772), 
for  which  the  aluminum  cylinder  had  a gauge  of  1.549  mm  (0.061  inches) 
and  the  lead  overlay  was  1.727  mm  (0.068  inches)  in  thickness.  As  shown 
in  Figure  2,  the  specimen  was  positioned  so  that  the  axis  of  the  electron 
beam  was  perpendicular  to  the  axis  of  the  cylinder  and  passed  through 
it.  In  future  references  to  circumferential  locations,  the  angle  6 will 
be  measured  from  the  direction  of  the  incident  electron  beam,  as  indicated. 
Foil  type  strain  gages  were  cemented  to  the  inside  of  the  cylinder  at 
angles  0°,  90°,  180°,  and  270°,  midway  between  the  ends  of  the  cylinder 


Figure  2.  Cylinder  Test  Configuration 
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with  their  axes  oriented  in  the  longitudinal  direction.  The  splice  in 
the  aluminum  cylinder  was  located  at  6 = 145°  while  the  splice  in  the 
lead  overlay  was  also  positioned  on  the  back  side  of  the  cylinder  within 
20°  from  0 = 180°. 

This  target  configuration  was  subjected  to  a pulsed  8 MeV  mono- 
energetic  electron  beam  with  a fluence  of  60  cal/cm2  in  the  drift  tube 
of  the  AURORA  facility.  Figure  3 is  a post-test  photograph  of  the 
aluminum  cylinder  which  shows  the  extent  of  permanent  deformation 
induced  in  this  specimen.  Quantitative  measurements  of  the  deformed 
cylinder  will  be  presented  when  comparisons  are  made  with  predictions. 
Figure  4 is  a photograph  of  the  recovered  soft  lead  overlay  in  a 
flattened  position.  On  the  average,  366  mm  (14.4  inches)  of  the  length 
of  this  strip  (corresponding  to  a circumferential  angle  of  212°)  was 
visibly  unaffected  by  the  electron  beam.  The  remaining  88°  of  the 
circumference  of  the  lead  cylinder  was  partially  melted  or  blown  off. 

At  each  end  of  the  residual  lead  overlay  there  is  a region  of  tapering 
thickness  approximately  64  mm  (2.5  inches)  wide  which  has  undergone 
metallurgical  transition.  Since  the  total  length  of  the  lead  strip  is 
appreciably  greater  than  the  initial  circumference  of  the  lead  overlay, 
the  affected  region  must  have  stretched  like  taffy  before  rupturing. 
There  is  no  evidence  that  the  electron  beam  impinged  directly  upon  the 
aluminum  cylinder. 


III.  MODELING  CONSIDERATIONS 

2 

The  previously  cited  study  of  the  coupling  of  x-ray  deposition  to 
structural  response  considered  the  problems  associated  with  linkage  of 

3 

two  analysis  models:  the  Radiation  Interaction  Program  (RIP)  and  the 

4 

PETROS  3.5  code  , a finite  difference  analysis  model  for  the  finite 
deflection,  elastoplastic  response  of  thin  shells  to  transient  loading. 
The  RIP  code  is  basically  a one-dimensional  hydrodynamic  stress  wave 
analysis  with  provision  for  introducing  various  radiation-related 
excitations.  For  x-ray  deposition  it  is  only  necessary  to  define  the 
source (s)  and  subroutines  perform  the  deposition  calculations.  However, 
for  electron  beam  excitation  the  provisions  are  much  more  rudimentary: 
it  is  necessary  to  specify  as  input  the  dose-depth  distribution  and  its 


' R.  H.  Fisher,  G.  A.  Lane,  and  R.  A.  Cecil,  "RIP,  A One-Dimensional 
Material  Response  Code, " Systems,  Science  and  Software  Reports 
3SR-751-I  <5  II,  September  1972. 

4 

S.  D.  Pirotin,  B.  A.  Berg,  and  E.  A.  Witmer,  "PETROS  3.5:  New 
Developments  and  Program  Manual  for  the  Finite-Difference  Calculation 
of  Large  Elastic-Plastic  Transient  Deformations  of  Multilayer 
Variable-Thickness  Shells,"  U.S.  Army  Ballistic  Research  Laboratories 
Contract  Report  No.  211,  February  1975.  (AD  #A007215) 
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time  dependence.  After  the  excitation  has  been  defined  the  RIP  code 
performs  finite  difference  calculations  of  stress  wave  propagation  and 
reflections  in  multilayer  plates.  When  the  solution  has  been  marched 
out  in  time  for  a sufficient  period,  information  may  be  extracted  from 
the  RIP  code  concerning  mean  velocities  acquired  by  the  structure, 
inelastic  stress  states  established  by  the  passage  of  intense  stress 
waves,  and  temperature  distributions  (which  may  be  sufficiently  high  that 
thermal  degradation  of  material  properties  must  be  considered)  for 
introduction  into  the  PETROS  3.5  code.  These  data  are  required  for  each 
mesh  point  of  the  two-dimensional  grid  defining  the  reference  surface 
of  the  shell  structure.  Where  this  input  information  varies  over  the 
surface  of  the  shell,  many  RIP  code  runs  may  be  required  for  an  adequate 
definition . 

The  authors  believe  that  the  weakest  link  in  this  modeling  chain 
is  between  the  output  of  the  field  emission  tube  and  the  input  to  the 
RIP  code  for  its  electron  beam  mode.  At  the  beginning  of  this 
investigation  the  experimenters  indicated  that  the  beam  is  perfectly 
uniform  with  an  average  angle  of  incidence  of  45°.  Since  prior 
experience  with  shell  response  codes  had  demonstrated  that  the  response 
pattern  of  the  cylinder  would  be  quite  sensitive  to  the  distribution  of 
surface  loading,  more  specific  deposition  information  was  requested.  We 
were  then  provided  with  predictions  of  dose-depth  distribution  obtained 

from  runs  with  the  one-dimensional  electron  transport  code  ZEBRA1"  for 
each  15°  of  circumferential  angle  of  the  front  side  of  the  cylinder 
(see  Figure  5). 

These  data  were  regarded  with  some  skepticism  since  there  was  no 
reduction  of  fluence  with  increasing  9 but  only  a shift  of  the  location 
of  maximum  energy  deposition  toward  the  outer  surface  of  the  lead. 
However,  these  data  were  employed  as  input  for  seven  RIP  runs  and  the 
results  were  interpolated  (where  necessary)  to  obtain  input  for  a 
PETROS  3.5  calculation  of  shell  deformation.  It  was  found  that  the 
RIP  runs  predicted  that  the  lead  was  either  melted  or  in  a solid- 
spalled  state  throughout  its  thickness  for  -90°  < 0 < 90°,  contrary  to 
the  evidence  of  Figure  4.  Also  the  PETROS  3.5  run  produced  a 
deformation  pattern  (shown  in  Figure  6)  which  differs  from  the  observed 
in  both  form  and  amplitude  (see  Figure  3 and  the  experimental  contours 
in  Figure  6) . It  must  be  concluded  that  we  had  not  employed  the 
correct  energy  distribution  as  input  for  the  RIP  code  runs. 

Based  partly  upon  the  observed  extent  of  damage  to  the  lead  overlay 
and  partly  upon  past  experience  in  predicting  cylinder  response  which 
indicated  that  a more  "peaked”  circumferential  variation  of  loading  was 
required  to  produce  the  pattern  of  Figure  3,  it  was  decided  to  employ 


L.  D.  Buxton , "The  Electron  Transport  Computer  Code  ZEBRA  1,"  Harry 
Diamond  Laboratories  TR-1536,  June  1971. 
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Figure  5.  Deposition  Profiles  for  Test  Configuration 


the  dose-depth  distribution  for  0=0  but  with  the  60  cal/cm'  fluence 
sealed  down  according  to  the  function 


J cos  20  for  I Ol  < 45° 

F(  0)  = 

( 0 for  45°  <r  I ol  < 180° 


f 1) 


It  will  be  seen  that  this  results  in  much  better  agreement  with  the 
available  test  data. 


IV.  THE0RETICA1, /EXPERIMENTAL  CORRELATIONS  FOR  THE 
IMPROVED  RESPONSE  MODEL 


A.  Material  Response 

The  dose-depth  distribution  resulting  from  use  of  Equation  (1)  is 
exhibited  in  Figure  7.  RIP  code  runs  using  the  electron  beam  mode  were 
made  for  each  of  the  angles  0 shown  in  this  figure.  Nine  degree  intervals 
were  chosen  to  match  the  nodal  point  locations  selected  for  the  PETROS 
3.5  code  calculation  to  follow.  For  each  RIP  code  calculation  the 
thickness  of  the  lead  overlay  was  subdivided  into  102  cells,  while  40 
cells  were  used  to  represent  the  thickness  of  the  aluminum  cylinder. 

The  GRAY  equation  of  state  was  employed  for  both  materials,  using  the 
properties  listed  in  Tables  I and  II. 

The  RIP  calculations  were  performed  in  the  normal  undamped  manner 
until  the  plastic  work  was  essentially  complete,  at  which  point  damping 
was  introduced  to  suppress  the  through-thickness  stress  waves.  This 
process  is  illustrated  in  Figure  8,  which  shows  the  variation  of  the 
kinetic  energy  of  a longitudinal  strip  of  aluminum  at  the  cited  angle 
(recall  that  this  is  a one-dimensional  calculation).  In  this  instance 
damping  was  initiated  slightly  after  1.4  ps.  The  residual  kinetic 
energy  is  that  of  uniform  radial  inward  motion.  The  corresponding 
variations  of  the  velocities  of  the  mass  centers  of  the  aluminum  strips 
at  various  angles  0 obtained  from  the  RIP  code  runs  are  depicted  in 
Figure  9.  When  damping  is  complete  there  is  no  through-thickness  motion 
relative  to  the  mass  center  and  the  stresses  in  this  direction  have 
vanished.  However,  when  plastic  work  has  been  done  there  will  generally 
exist  residual  tangential  stresses;  i.e.,  in-plane  stresses  in  the 
cylindrical  shell.  Figure  10  shows  a representative  through-thickness 
distribution  of  these  residual  tangential  stresses,  which  must  be 
transferred  as  input  to  the  structural  response  calculation.  Also 
obtained  from  the  RIP  code  solutions  was  information  on  temperature 
distributions  produced  by  the  electron  beam  deposition,  a typical  one 
of  which  is  shown  in  Figure  11.  Since  for  all  angles  the  aluminum 
experienced  a temperature  rise  of  only  a few  degrees  above  ambient 
(300°K) , it  may  be  concluded  that  it  is  unnecessary  to  give  further 
consideration  to  thermal  stresses  or  thermal  degradation  of  material 
properties  in  the  response  modeling. 
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Figure  7.  Modified  Deposition  Pro 


Table 

I . 

ciUAY  Equation  of  State  Parameters  for  Lead 

Symbo 1 

Numerical  Value 

Property 

po 

1 1.355 

gm/cm2 

Ambient  density 

t: 

0 

2.1  x 

105 

cm/sec 

Bulk  sound  speed 

AMU 

5.56  x 

10l° 

dynes/cm2 

Shear  Modulus 

Y 

0 

1 .0  x 

108 

dynes/cm2 

Yield  strength 

0 

o 

-1.0  x 

10l° 

dynes/cm2 

Spall  limit 

l:vv 

9.43  x 

109 

ergs/cm 

Material  has  vaporized 

f 

LV 

2.90  x 

109 

ergs/cm 

Material  commences  to  vaporize 

E 

LM 

X 

o 

C'j 

xO 

108 

ergs/cm 

Material  has  completed  melting 

F 

SM 

4.10  X 

10r 

ergs/cm 

Material  commences  to  melt 

S 

1 .54 

) 

Hugoniot  parameter 

Yo 

2.84 

dimensionless 

Lattice  gamma 

a 

2.3 

\ 

a = y - . 5 

0 

Ye 

. 6667 

' 

Electronic  gamma 

gc 

14.7  x 

10-9 

Mbar-cm3/mole-deg2 

Electronic  energy  coefficient 

*mo 

760. 

K 

Melting  temperature  parameter 

AW 

207.2 

gm/mole 

Atomic  weight 

vj 

vb 

.109 

cm3/gm 

Volume  at  which  EOS  are  joined 

.042 

cm3/gm 

Excluded  volume  for  vapor  phas 

a 

49. 

Mbar  (cm3/mole) 

Coefficient  of  attractive 

y 

potential  for  vapor 

e 

1.0 

dimensionless 

Joint  parameter 

DELS 

9.637 

X 10- 

5 Mbar-cm3/mole-deg 

Entropy  of  melting 

NOTH:  Cgs  rather  than  SI  units  are  listed  in  Tables  I and  II  since  these 
units  are  required  for  the  present  version  of  the  RIP  code. 
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Table  II.  GUAY  Equation  of  State  Parameters  for  6061-T6  Aluminum 


Symbol 


Numerical  Value 


po 

2.704 

gm/cm3 

Ambient  density 

Co 

5.21  x IQ5 

cm /sec 

Bulk  sound  speed 

AMU 

2.79  x 1011 

dynes/cm2 

Shear  Modulus 

Y 

0 

2.77  x 109 

dynes/cm2 

Yield  strength 

a 

0 

-1.3  x 1Q10 

dynes/cm2 

Spall  limit 

Evv 

1.2  x 1011 

ergs/cm 

Material  has  vaporized 

elv 

3.05  x 101 

0 

ergs/cm 

Material  commences  to  vaporize 

-u. 

1.1  x 1010 

ergs/cm 

Material  has  completed  melting 

E 

SM 

7.15  x 10C 

ergs/cm 

Material  commences  to  melt 

S 

1.338 

Hugoniot  parameter 

Y 

2.18 

Lattice  gamma 

'o 

f 

dimensionless 

a 

1 .68  | 

a = y - . 5 

0 

.6667  i 

Electronic  gamma 

ge 

8.7  x if)-9 

Mbar-cm3/mole-deg2 

Electronic  energy  coefficient 

T 

mo 

1340. 

K 

Melting  temperature  parameter 

AW 

26.98 

gm/mole 

Atomic  weight 

VJ 

.452 

cm3/gm 

Volume  at  which  EOS  are  joined 

Vb 

. 190 

cm3/gm 
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Figure  8.  Kinetic  Energy  of  Aluminum  Strip 
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DEPTH  MEASURED  FROM  INITIAL  OUTER  SURFACE,  (mm) 


The*  RIP  code  solutions  also  yield  information  concerning  the 
thickness  of  solid  lead  remaining  after  electron  beam  deposition  as 
a function  of  circumferential  angle  (see  Figure  12).  This  variation 
of  thickness  appears  to  be  in  reasonable  agreement  with  that  of  the 
recovered  overlay.  Precise  thickness  comparisons  are  not  feasible  since 
the  overlay  experienced  stretching  and  additional  thinning  before 
rupturing.  Also  plotted  in  Figure  12  (to  an  arbitrary  scale)  is  the 
circumferential  variation  of  radial  velocity  of  the  aluminum  cylinder 
which  is  to  be  transferred  as  an  initial  condition  to  the  structural 
response  analysis. 


Figure  12.  Thickness  and  Velocity  Distributions 
B.  Caveats  Concerning  Imperfect  Modeling 

Before  proceeding  with  consideration  of  response  predictions  it  is 
appropriate  to  note  those  non-ideal  aspects  of  the  experimental 
procedure  which  will  be  ignored  in  the  modeling  but  which  may  contribute 
to  discrepancies  between  predictions  and  experimental  data.  By 
examination  of  the  recovered  lead  overlay  (see  Figure  4)  it  may  be 
inferred  that  the  electron  beam  deposition  was  not  completely  uniform 
for  the  longitudinal  direction  of  the  cylinder  and  that  the  variation 
of  intensity  was  not  symmetric  with  respect  to  the  mid-length  of  the 
cylinder.  This  effect  is  apparent  in  the  deformed  cylinder  (Figure  3) 
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also,  whore  the  maximum  deflection  is  located  about  12.7  mm  (0.5  inch] 
from  the  mid-cross  section.  Further,  it  is  observed  that  the  response 
pattern  of  the  cylinder  seems  to  be  rotated  about  12°  with  respect  to 
the  9 = 0°  plane.  There  are  several  possible  causes  for  this  effect: 

(a)  The  9 = 0°  strain  gage  may  not  have  been  aligned  with  the 
incident  beam. 

(b)  The  incident  beam  axis  may  have  been  above  or  below  the 
cylinder  axis. 

(c)  The  looseness  of  the  fit  of  the  cylinder  in  the  end  caps  may 
have  permitted  a rotation  when  the  cylinder  was  struck  by  the  beam 
(pendulum  effect) . 

(d)  Geometric  or  material  non-uniforunities  of  the  specimen. 

It  does  not  appear  possible  at  this  time  to  pin-point  one  of  these 
possible  causes  as  being  responsible.  Another  effect  associated  with 
the  looseness  of  the  end  caps  is  the  "pulling-in"  of  material  in  the 
longitudinal  plane  of  maximum  deflection  and  the  flaring  out  of  material 
which  overlays  the  end  caps.  As  a measure  of  this  effect,  the  total 
length  of  the  cylinder  prior  to  deformation  was  190.5  mm  (7.5  inches), 
whereas  after  deformation  the  length  at  the  region  of  maximum  transverse 
deflection  was  reduced  to  185.4  mm  (7.3  inches).  The  effect  of  this 
"pulling-in"  of  material  is  to  permit  larger  transverse  deflections  than 
would  occur  if  this  were  inhibited.  A converse  effect  resulted  from 
an  error  in  communication : after  all  calculations  were  completed  it  was 
discovered  that  the  aluminum  cylinder  was  actually  0.051  mm  (0.002 
inches)  thicker  than  the  value  previously  cited.  Thus,  the  cylinder 
should  be  somewhat  stiffer  than  predicted,  especially  with  regard  to 
flexure. 

C.  Structural  Response 

The  structural  response  prediction  was  performed  using  the  PETROS 
3.5  code,  the  aluminum  cylinder  being  modeled  as  having  fixed  ends  and 
as  being  fabricated  from  an  elastic,  perfectly  plastic  material  (see 
Figure  1)  with  a yield  stress  of  276.5  MPa  (40100  psi)  and  a density  of 
2701  kg/m3(0. 0002528  lb  sec2/in4) . The  lead  overlay  was  ignored  in 
this  calculation  since  (a)  it  was  not  bonded  to  the  aluminum  cylinder 
and  (b)  high  speed  photography  showed  the  overlay  rapidly  peeling  away 
from  the  structural  shell.  The  distributions  of  residual  tangential 
stresses  and  the  velocity  distribution  of  Figure  12  were  used  as 
initial  conditions  for  the  response  calculation,  assuming  that  there 
was  no  longitudinal  variation  of  these  quantities. 

The  evolution  of  the  final  deformation  pattern  is  illustrated  in 
Figure  13.  The  grid  defining  the  cylinder  surface  in  the  isometric  views 
is  the  actual  finite  difference  grid  employed  for  the  calculations.  The 
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Figure  13.  Transient  Deformation  Patterns  of  the  Aluminum  Cylinder 


deformation  of  the  cylinder  is  plotted  to  the  same  scale  as  the  initial 
geometry;  i.e.,  the  deflections  are  not  exaggerated.  A larger  view  of 
the  predicted  final  deformation  mode  is  shown  in  Figure  14.  For 
comparison,  the  post-test  deformation  pattern  measured  from  the 
cylinder  of  Figure  3 is  also  plotted  in  the  sectional  views. 

Considering  the  rather  arbitrarily  adapted  load  distribution  and  the 
numerous  caveats  necessary  with  respect  to  imperfections  in  the 
modeling,  the  agreement  between  theory  and  experiment  is  regarded  as 
quite  satisfactory. 

The  predicted  transient  displacement  of  the  point  on  the  crown 
line  midway  between  the  ends  of  the  cylinder  (labeled  point  B on 
Figure  14)  is  displayed  in  Figure  15.  It  is  seen  that  the  deflection 
at  this  position  of  maximum  deformation  increases  monotonically  until 
the  maximum  is  reached,  with  little  recovery  ensuing.  This  is 
consistent  with  the  energy  balance  diagram  of  Figure  16,  where  it  may  be 
seen  that  the  plastic  work  is  essentially  complete  by  260  as,  the 
subsequent  small  amplitude  oscillations  being  mainly  elastic.  The 
small  gap  between  the  upper  two  curves  of  this  figure  are  an  indication 
of  the  discretization  error  of  the  numerical  analysis. 

The  predicted  transient  strain  components  at  each  90°  position 
(Points  B,  A 5 C,  D of  Figure  14)  are  shown  in  Figures  17,  18,  and  19, 
respectively.  One  sees  the  expected  delayed  arrival  of  strains  at 
6 = 90°  and  180°.  Also  the  strains  are  quite  large  for  9=0°,  with 
membrane  strains  dominant  over  bending  strains.  At  9 = ± 90°  the 
strains  are  much  smaller  (but  entailing  some  plasticity),  beginning  as 
principally  membrane  strains  and  then  undergoing  a transition  to  mostly 
bending  strains.  The  strains  at  0 = 180°  are  of  the  same  order  of 
magnitude  as  those  at  6 = 90°  but  undergo  less  of  a transition  to 
bending  strains.  These  predicted  strains  are  entirely  consistent  with 
the  deformation  patterns  of  Figure  13. 

It  had  been  hoped  to  display  the  experimental  data  for  the 
longitudinal  strains  on  the  inner  surface  of  the  cylinder  on  Figures 
17,  18,  and  19  for  comparison  with  the  model  predictions.  However,  by 
reference  to  Figure  20,  which  includes  reproductions  of  the  oscilloscope 
pictures  taken  during  the  test,  one  sees  that  this  is  hardly  feasible. 

It  may  be  seen  that  the  horizontal  sweep  speeds  employed  were  too  slow 
to  permit  observations  of  details  of  the  strain  fluctuations  provided 
by  the  PETROS  3.5  prediction.  The  alternative  possibility,  extending 
the  PETROS  3.5  solution  for  the  duration  of  the  oscilloscope  records 
is  not  desirable,  both  due  to  the  excessive  computer  time  required  and 
because  no  attempt  has  been  made  to  model  the  internal  damping  in  the 
structural  cylinder.  Another  difficulty  with  the  experimental  data 
concerns  the  location  of  zero  for  the  vertical  axis;  this  was  initially 
set  on  the  horizontal  axis  having  tic  marks  but  may  have  drifted 
somewhat  before  the  electron  beam  pulse  was  activated.  Further,  no 
trace  was  recorded  for  gage  A and  both  gages  A and  B were  found  to  be 
disconnected  after  the  test.  It  is  believed  that  gage  B failed  very 
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Figure  16.  Energy  Balance  Diagram  for  Aluminum  Cylinder  Response 
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Figure  1".  Strains  at  0 0°  (Point  B) 
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Figure  19.  Strains  at  0 = 180°  (Point  P) 
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early  during  the  test  since  the  recorded  trace  amplitude  is 
insignificant  in  comparison  to  the  predicted  maximum  strain  of 
approximately  9%. 

For  gage  C the  correlation  between  theory  and  experiment  is  rather 
good:  the  predicted  maximum  strain  is  0.2%  while  the  experimental  value 
is  at  least  0.17%,  perhaps  more  if  the  trace  went  off  scale.  Unfortu- 
nately, the  experimental  trace  for  gage  D is  invisible  against  the 
dense  background  for  the  first  seven  milliseconds,  during  which  the 
largest  strains  are  predicted. 

In  an  effort  to  predict  the  residual  strains  in  the  aluminum 
cylinder  the  PETROS  3.5  calculations  were  restarted  with  damping  and 
run  until  the  kinetic  energy  became  negligibly  small.  These  predicted 
residual  strains  are  compared  to  the  corresponding  experimental  values 
in  Table  III.  The  agreement  between  prediction  and  test  for  gages  C and 

Table  III.  Residual  Longitudinal  Strains  in  Aluminum  Cylinder 


0 Gage 


Longitudinal  Strain 
Predicted  Experimental 


0°  B 

90°  C 


8.29%  invalid 

0.04%  0.01% 


180°  D 


-0.11%  -0.09% 


D may  be  seen  to  be  quite  good  when  it  is  appreciated  *hat  the 
discrepancy  is  small  in  comparison  to  the  magnitudes  of  the  transient 
strains  at  these  locations.  While  it  is  dissappointing  that  a 
graphical  comparison  of  predicted  and  measured  transient  strains  is  not 
available,  it  can  be  stated  that  no  experimental  strain  data  which 
appear  valid  are  in  significant  disagreement  with  predictions. 


V.  CONCLUDING  REMARKS 

It  is  believed  that  the  objective  of  demonstrating  the 
capability  of  a combined  ZEBRA/RIP/PETROS  3.5  analysis  to  model  the 
response  of  a generic  shell  structure  to  electron  beam  loading  has 
been  achieved.  The  correspondence  between  predictions  and  test  data 
was  adequate  for  the  quality  of  experimental  results  available.  In 
future  tests  even  better  correspondence  should  be  attainable  by 
devoting  greater  care  to  control  or  measurement  of: 

(a)  gage  thickness  and  geometrical  imperfections  of  specimens 

(b)  material  properties  (including  possible  rate  effects) 

(c)  boundary  conditions  (which  can  be  readily  modeled). 
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The  performance  of  pre-test  predictions  should  assist  in  the  selection 
and  adjustment  of  instrumentation  and  result  in  better  quality 
experimental  data. 


As  noted  previously,  the  only  obstacle  to  successful  use  of  the 
AURORA  facility  in  the  electron  beam  mode  as  an  x-ray  loading  simulator 
is  the  apparent  absence  of  an  analytical  model  relating  the  electron 
beam  "free  field"  to  the  fluence  distribution  over  the  surface  of  the 
target;  i.e.,  to  the  input  to  the  ZEBRA  calculations.  This  modeling  gap 
was  circumvented  in  the  present  investigation  by  an  intuitively  deduced 
circumferential  distribution  which  must  have  been  reasonably  close  to 
the  physical  distribution,  judged  by  the  success  in  predicting  the 
deformation  pattern  of  the  target.  However,  for  future  more  complex 
applications  a satisfactory  model  for  fluence  distribution  should  be 
developed . 
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